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Explicit and implicit schemes for 
fractional—order Hantavirus model 
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Abstract 


In this paper, the fractional-order form of a mouse population model is 
introduced. Some explicit and implicit schemes, which are Theta methods 
and nonstandard finite difference (NSFD) schemes, are implemented to give a 
numerical solution of nonlinear ordinary differential equation system named 
Hantavirus epidemic model. These methods are compared and discussed that 
the method preserves the positivity properties of the integer order system. 
The numerical solutions are illustrated by means of some graphs. Numerical 
results of explicit and implicit methods denote that these methods are easy 
and accurate when applied to fractional—order Hantavirus model. 


Keywords: Explicit and implicit methods; Theta method; Nonstandard fi- 
nite difference scheme; Fractional-order nonlinear differential equation sys- 
tems; Mouse population model. 


1 Introduction 


Although fractional calculus has a long history, its applications to natural sci- 
ence are just a recent focus of interest. The description of some phenomena 
is more accurate when the fractional derivative is used. In many scientific 
area specially in physics, chemistry, and engineering, the fractional differ- 
ential equations (FDEs) become a popular subject, which are increasingly 
used to model problems in a number of research areas including dynamical 
systems, mechanical systems, signal processing, electronic circuit theory, con- 
trol theory, chaos synchronization, mechanics, seismology, and many others. 


*Corresponding author 

Received 19 September 2015; revised 18 April 2018; accepted 16 May 2018 
Mevliide Yakit Ongun 

Department of Mathematics, Suleyman Demirel University, Isparta, Turkey. 
e-mail: mevludeyakit@sdu.edu.tr 


Damla Arslan 
Graduate School of Natural and Applied Sciences, Suleyman Demirel University, Isparta, 
Turkey. 


75 


76 Mevlide Yakit Ongun and Damla Arslan 


Some of these studies may be found in [4] and [B5). Lots of books written 
by the authors Podlubny [80], Lubich [19], Miller and Ross [24], Oldham and 
Spainer [28], Diethelm et al. [8], Samko et al. [B4] played a significant role to 
understand the subject of FDEs and gave some methods for solutions. 

Several methods, which are the Adams—Moulton method [14], [2], the 
homotopy perturbation method [25], the Adomian decomposition method 
[25], [26], the variational iteration method [B6], the differential transform 
method [8], [Ld], the operational method [20], the predictor corrector meth- 
ods [@, (8), [3], the product integration rules [15], [B83], nonstandard finite 
difference scheme [29] have been used to solve FDEs. Numerical solution 
methods are frequently referred to as being explicit or implicit. These sit- 
uations are being encountered numerical solutions of FDEs and equation 
systems. 

The first and important study, which is relevant with Hantavirus model, 
is given by [2]. The stability analysis of model, which is discretized with 
nonstandard finite difference scheme, is given by [5], [6]. Further, variational 
iteration method [16], exponential matrix method (EMM) [89], multistage 
differential transformation method [IZ] have been used to solve Hantavirus 
model and FDE solutions of this model studied by [81], [Bz]. 

The effect of the environmental parameter on the Hantavirus infection 
through a fractional-order SI model was studied by Rida et al.(2012). They 
presented a fractional-order model of the Hantavirus infection in terms of 
simple differential equations involving the mice population. They studied 
that the effect of changes in ecological conditions and diversity of habitats 
can be observed by varying the value of the environmental parameter k. They 
used a generalized Euler method (GEM) to obtain an analytic approximate 
solution of the model [87]. 

In this paper, new numerical methods determined for fractional—order 
nonlinear differential equation model known as Hantavirus epidemic model. 
These model based on a nonlinear differential equation of order p, where p is 
a constant in range 0 < p < 1. Some explicit and implicit methods such as 
theta method and nonstandard finite difference schemes are studied for the 
numerical solution of the fractional-order model. Especially the purpose of 
the studying to these two methods that both methods are more general face 
of classical methods and both methods allow to some arbitrary choice. 

The following model is given in [2] has been used in the study of Han- 
tavirus epidemics: 


dM, M,M 
p= OM ~ cM, aM.My, 
dM; M,M 
=~0M. M,Mz, 1 
dt i aie oR a ” 


where 
M, : the population of susceptible mice (M, > 0), 
M, : the population of infected mice (My; > 0), 
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b: the birth rate, 

c: the death rate, 

K : the carrying capacity of the environment. 

In this model, mice movement as a process of diffusion is ignored and 
whole population is composed of two cases of mice, susceptible and infected. 


For the systems (ll), the total population M = M,+ M, satisfying the logistic 
differential equation 


=(b-M-—. (2) 


The carrying capacity for the total population is M* = (b—c)K. The 
critical value of the carrying capacity is 


Ke= 


a(b—c) 


Let us M, = x, My = y in equation (1), 


dy _ y? l-ak 3) 
de K me 


The dynamics of the continuous system (B) have been given in [2,,5,,6] as 
the following: 

i) If b < c, then the system has a unique equilibrium (0,0), and it is 
globally asymptotically stable. 

ii) fb >cand K < COE then the system has two equilibria: (0,0), 
which is unstable, and (A’(b—c),0), which is globally asymptotically stable. 

iii) Ifb > cand Kk > ACE then the system has three equilibria: (0,0) 
and (K(b—c), 0), which are unstable; and (2, K(b-—c)- b) , which is globally 
asymptotically stable. 

In this article,the fractional-order Hantavirus model examined, which 
given as below: 
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where noninteger derivative is defined as the Caputo derivative. The fractional— 
order equation system is introduced and denoted that this model has a unique 
solution with given initial conditions while t > 0. Some explicit and implicit 
schemes are implemented to study the dynamics behaviors of fractional order 
Hantavirus epidemic system. Some numerical solutions illustrated by means 
of some graphs are provided in next sections. 


2 Fractional derivatives 


On this section, some definitions and relationship of fractional order integra- 
tion and fractional-order differentiation mentioned [29]. 


The Griinwald—Letnikov (GL) operator of order p > 0 is defined as 


N 
ouDi, F(t) = fim hy? Dw)" F(t Shy) 
j=0 


w!?) = (-1)! @ __ TG») 


with 
j Diet +i) 


Note that the weights wi? ) are the coefficient in the power series expansion 
of (1 — €)?; that is, 


co 
-9r= lure 
j=0 
and they can be evaluated recursively by means of the recurrence 


1 
tea, w= (1-242) ul, FH... (4) 


The relationship between Riemann-Liouville (RL) and Caputo definitions 
is 


oD, f(t) =rz Di, (f(t) — Tm-iIF; tol) ; 


where Tj,—1[f;to] is the (m — 1)th degree Taylor polynomial for f centered 
at to 
m—1 1k 
Tm—1 [f; tol (¢) — an (to). 


k=0 


While, 0 < p < 1, it ism =1 and To[f; to|(t) = f(to). Thus, we obtained 


oD? f(t) =ri D?, (f(t) — f(to)). 
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Under suitable assumptions of regularity, the RL and GL operators coin- 
cide; that is, 


riDi, f(t) =er Dif). 


A consequence of the above two points is that in general 


cDi, f(t) =ax Di, f(O) (F(t) — Tm-alf; to]) ; 


and while 0 < p < 1 it is 


oD, f(t) =er Do.(F(t) — f(to))- (5) 


In case it used in differential equations of fractional order, the Caputo 
definition is preferable since it allows the couple the equation with initial 
conditions of classical type (i.e., initial condition of Cauchy type). In this 


way we can obtain 
y(to) = Yo- 


In case we use the RL definition; we can couple different initial conditions, 
as 


lim Ju, y(t) = bj, 


+ 
totd 


{ini = f(t,y@), 


which does not have a clear physical meaning. Therefore, they are not 
useful for practical applications. 


Lemma 1. Let 0 < p< 1, and let wl?) be the coefficients in the power series 


expansion of (1 — €)% given in w?). Then for anyn =1,2,..., 


J 
a) -1<w”) <0, 


bO<w?) <1. 


Proof. The proof is an immediate consequence of the recursive relationship 
stated in (@). 


3 Fractional order Hantavirus epidemics 


In this section, the fractional order applied in the model of Hantavirus epi- 
demics [2]. The new system is given as the set of fractional differential equa- 
tions (FDEs): 
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Or (Fe) cout, 


eX) tuto, (6) 


«(0)=29, y(0)= yo, (7) 


where cD}, denotes the fractional derivative operator, with respect to the 
origin, according to Caputo’s definition [BD] and x(t) and y(t) are activa- 
tor and inhibitor variables. The fractional derivatives are used to describe 
nonhomogeneous character of the ecosystems, with respect to the presence of 
competitors. The parameter p denotes the density of competitor species in the 
system. The reason, to use fractional order differential equations (FODE), 
is to be naturally related with systems with memory, which exists in most 
biological systems. Also FODE are closely related to fractals, which are abun- 
dant in biological systems [7]. The results derived of the fractional system 
(B) and (@) are more general nature. 

When the power exponent is p = 1, this corresponds to equation (8) and 
varies competitor’s populations when 0 < p < 1. While p > 1, the density of 
competitor or alien species will increase in the populations [I]. The stability 
analysis given by the study of [29]. 

For non-negative solutions, denote RZ = {x €R?|a>0}, and let 
a(t) = (X,Y)". To prove the main theorem, we need following lemma and 
corollary given by [9,27]. 


Lemma 2 (Generalized Mean Value Theorem). Suppose that f(x) € Cla, }] 
and D? f(a) € Cla, b] for 0 < p< 1; then we have 


1 P —q)P 
f(x) = f(a) + T® (Dif) (6)(@ — @) 


witha <€ <a for all x € (a,b). 


Corollary 1. Suppose that f(x) € Cla,b] and D? f(x) € Cla,b] for 0 < 
p<. If D’ f(x) = 0, for all x € (a,b), then f(x) is nondecreasing for each 
x € [a,b]. If DP f(x) < 0, for all x € (a,b), then f(x) ts nonincreasing for 
each x € [a,b]. 


Theorem 1. There is a unique solution x(t) = (X,Y)! to equation (8) 
with initial condition (Q) on t > 0, and the solution will remain in R%.. 
Furthermore, X and Y are all bounded by Kg. 


Proof. We know that the solution on (0,00) solving the initial value problem 
(G)-(@) is not only existence but also unique [I8}. 
From equation (@), we find 
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D?X |x=0= bY > 0, DY |y=o= 0. 


By Corollary 1, the solution will remain in R%. From [2], X and Y are 
bounded by Ke. 


To evaluate the equilibrium points, let DPx = 0 and D?y = 0. Then 
Eo = (0,0), Ei = (K(b—c),0), andE2 = (2, K(b—c)— 2) are the equilibrium 
points. To give more detail about the local behavior near the equilibria, we 
find the Jacobian matrix of equation (@) at each equilibrium point: 


B= (7° 2) sEd= (Ne? tea) 


and 


_ (-aK(b-¢)+0(1- 4) — 
HE) = (cee noe ee 


a 


The equilibrium point Eo is locally asymptotically stable, if all of the 
eigenvalues A; (i = 1,2) of the Jacobian matrix J(F;) for i = 0,1,2 satisfy 
the following condition [10,21]: 


larg(As)| > (8) 

The eigenvalues of J(Fo) are Ay = (b—c) and Az = —c. If b > c, then 
arg(A,) = 0, arg(A2) = 7, and Ep is unstable. Ifb < c, arg(\,) = arg(A2) = 7, 
and equation (§) is correct, then Eo is asymptotically stable. 

Now consider the case b > c and the equilibrium point FE. The eigenvalues 
of J(£) are Ay = —(b—c) and Ap = —b+ ak (b—c). arg(A1) = 7 and for 
K < Ke, arg(A2) = 7, thus E) is asymptotically stable. If kK > K., then Fy 
is unstable. 

Finally, consider the case b > c for the equilibrium point Fo. If K > Ke, 
then Ay = —(b—c) and Ay = b — ak (b —) are negative eigenvalues and 
arg(A1) = arg(A2) = 7; consequently, Ey is asymptotically stable. 


4 Theta method for fractional order models 


In this section, theta method (6 — —method) is described and studied for 
the numerical solution of fractional-order equation systems by the way of 
Hantavirus model. Theta method is known as the weighted method and 
general form given by 


DP? y(tn) = OF (tr, Yn) + (1 — 8) fF (tn—1; Yn—1) (9) 
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while @ € [0,1]. There are two specific values of 6: while 6 = 0, we obtain 
the forward explicit Euler method and while 6 = 1 yields forward implicit 
Euler method. Note that @ = $ in (@) corresponds the so-called mid-point 
rule and trapezoidal rule, respectively. 

For numerical computation the Griinwald—Letnikov (GL) operator is trun- 
cated; we fixed a step-size h > 0 and t, = to thn, N = (T — to)/h. 

We put the left-hand-side with the forward discrete derivative with 


aCe _ yo) = Of (tr, Yn) 1 (1 3 0) f (tr—15 Yn—1)- 
j=0 


Case 1. While 6 = 0 on (Q), in case we use it for the model equation of 
(8), we can get 


nr 2 
h-P Ye w§ (2n—j — Xo) = (b— €)&n—1 + bYn-1 a (+424) Ln—-1Yn-1; 
j=0 
n 2 
h-? Y) wh (Yn—7 — Yo) = —CYn—1 — “5g? — (4752*) an-1)-1 


Here, the left-hand derivatives are Griinwald—Letnikov derivatives. 


nm p n p 
S wlitnj — m0 oul 
j=0 j=0 
2 
7 


= hP le 2) teat bye a — (44) Zn—19n-1| , 


n n 2 
Yo whan 3 — wo Lowe = AP | cna — Yet — (159K) on ty] - 
j=0 = 


a2 
igi lo C)En—1 + bYyn—1 — 34 (42%) on1ym—a| 
-»> we? an F tw ao, 
j=l 

i —a 7 —1 

Yn = hP CYn—-1 fant (+ a) In 1Yn-1| _ 2 we? yn —j _ {p yo 
jJ= 
(10) 

Since wh” = —p, we obtain 
In-1 = we Yao, n=1 and wP Va, oe we? tn, n> 2, 


Explicit and implicit schemes for fractional-order Hantavirus model 83 


and similarly we can say for Yn—1. While 


(p—1) 


we ine — Xn =Fn-1 + PLn-1, (11) 
wl! Yo — Y, = Yn—1 + PYn-1; 


and 


We get (M1) in the equations system (10), and we obtain iteration equa- 
tions as: 


In = Ln—1 + Pln—1 
2 
+h? [@- C)@n—1 + byn—1 — 4 — (48 ) 2n—19%-1| > (12) 


Yn = Un-1 + PYn-1 + hP CYn-1 ue) (Hk 


Case 2. While 9 = 1 on (Q), use it for Equation (@), we obtain implicit 
form 


n 
h-P » wy (fn—j x0) = (b — C)in a byn ~ a _ (+424) TnYn> 
j=0 
mt 2 
h-P » wi (Yn—5 Yo) = —CYn i (43) TnYn- 


j=0 


In case we leave x, and y, alone on the left side, 


Ln hP (b C)Ln t byn ua (+424) nn] + Xn _ wl) x9 = 0, 


Yn — bP |—cyn, 7 (42*) nn] + Yn — we yo = 0. 


While system has implicit form, Newton—Raphson method needed to solve 
implicit system by converting to explicit form. Here, the implicit form, x, 
and yp, connected to the explicit form of the first equation to be able to resort 
to Newton—Raphson method: 


f(z) =24 Xp, —we Ya — bP [(e- c)z + bu — e — (+3*) zu ; 
gu) =utYn — we yo hP cu ue (454) zu ; 
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Regarding to Newton —Raphson method, iterative equations as: 


f(z) 

4 = ’ ’ 1, seg h Gy 
ates e f(z) : ° 
Uit1 = Ui — glu) ’ 0, 1, vey 

g' (ui) 


Here, we can write initial estimations as z9 = @,_1 and ug = Yn_1. Thus 
we obtain 


24—-En—1—Pz0 ne[(o c)zitbui- 


ee 1—he [(b—)— 2 —( FJ] 
2 (13) 
Ui—Yn—1—puo—h? | —cui— zB (45 )ziu: 
test = Tee EE] 


5 Nonstandard difference schemes methods for 
fractional order Hanta model 


In this section, we proposed and discussed some NSFD schemes to discretize 
the fractional order nonlinear system (@). NSFD schemes applied in combi- 
nation with the truncation of the GL operator. NSFD schemes were firstly 
proposed by Mickens for Ordinary Differential Equations and successively 
studied, for instance, in [22,23]. The stability analysis of model which is 
discretized with nonstandard finite difference scheme is given by [5,|6]. 

Case 1: As a first nonstandard scheme, for discretization we make the 
replacement of the nonlinear term in the right-hand side of (@) by means 
of a(t) + 2(tr—1), y(t) + yltn—1), a2(t) + x(tn)a(tn1), 2(t)y(t) > 
t(tn)y(tn—1)- 

Since w} = 1, the discretized model is 


In + Xn — we tx = o(h)X, 
Yn + Yn . we-tyo = P(h)Y, 


which seem as explicit form and ¢(h) is a denominator function and could 
* 


be selected in a different way. Here, X,,,Y,,X and Y are defined as: 


Xp = Sw? ens, Y, = Sw yn is 
j=l j=l 
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(@n—1+Yn-1) Yn 


. (he b (@n—-1+Yn-1)@n 
X = (b— c)&p—-1 + byn-1 ik ALnYn—-1; 
OK 
Y= CYn-1 K alnYn-1; 


_ K($(h)(b-c)n—1+(h)byn—1—Xn +P te) 


In = K+@(h)(@n—1+¥n—1)+K O(h)ayn—1 : 
= K(—$(h)cyn—-1—44(h)a@n¥n—1—-Yn tw? yo) 
a K+6(h)(@n—14Yn—-1) 


The result of explicit nonstandard finite differences scheme is 


K (@n-1+(p—cd(h))an—1+bG(h)(an—-1+Yn-1)) 
K+6(h)(@n-1+yn—-1)+aK $(h)yn—1 ; 


Ln = 


(14) 
— KGn-1t(p—¢h(h)) yn-1+40(h)@nYn-1) 
Yn = K+6(h)(@n—-1+Yn—1) , 


p 
where we selected ¢(h) = (aes) . In [82,83], authors have chosen 


—e 


(xy) + x(t+h)y(t), 2? > a(t)x(t+h) and y? > y(t)y(t+h) for discretization. 


Proposition 1. Let a,b,c, K > 0, b<c, and let step size 0 < h < (p/c)“/”) 
for the schemes have the iterations 2, Yn, given by the non-negative schemes 
for any xp > 0 and yo = 0. 


Proof. We proceed by induction on n. 


Case 2: In the second nonstandard scheme we use the replacement 
a(t) + (tn), y(t) > y(tr), e°(t) 4 2(tn)e(tr), a(t)y(t) + c(tr)y(tn). 
By operating in a similar way as in the previous case for system () , we can 
see that is implicit form: 


4k 


In + Xn —we-t29 = o(h)X, 
Yn + Yn = we lyo = P(h)Y, 


X = (b—c)an + dyn — Keates — a2nYn; 
(an+Yyn)Yyn 
K 


= —CYn aALnYn- 


If we leave alone to x, on the left side and y,, on the right side, we obtain 
implicit form: 

Ln — (h)X + Xn — we Ya = 0, 

Un — (h)Y +¥q —w? Yyp =0. 
While Newton—Raphson method is needed, here, the implicit form, x, and 


Yn connected to the explicit form of the first equation to be able to resort to 
the method of Newton—Raphson, we leave alone on the left side: 
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f(z) =2+Xn- wPYyp — o(h) |(b—c)z + bu — a - (+424) zu ; 
g(u) = Ut Yn — wi? up = b(h) [cu = 42 = (GK) zu]. 


Regarding to Newton—Raphson method, the initial estimations as zo = 
Ln—1 and up = Yn—1. In this way, we can obtain 


2 
2i—Bn—1—pz0—P(h) | (b—c) zi +bui z (+28) anu 


1—4(h)[(b—0)— 4 — (Eu ] 


Al = % 


are 
Ui -Yn-1—PUo (n)| cui— (52) 


1—$(h)[—e— Ht — (258% )ai] 


Ui41 = Ui 


6 Numerical simulations 


In this section, numerical simulations of fractional-order Hantavirus model 
are presented by the solution of Theta method and NSFD schemes that are 
explicit and implicit form. We choose a = 0.1, b= 1, and c= 0.5 as used by 
study of [2]. The initial values is considered as (ao, yo) = (25,15); we fixed 
a step-size h = 0.01 and N = 2000. Various denominator functions listed 
in [29], we took into account the ¢(h) = expt)" We choose p = 1 
in Figure I, p = 0.6 in Figure 2, p = 0.3 in Figure B p = 1 in Figure @, and 
(to, Yo) = (10,10) and K = 20 used in Figure &. In Figures [4], the top 
graphs in order of left-to-right are theta explicit and theta implicit methods, 
the bottom graphs in order of left-to-right are NSFD explicit and NSFD 
implicit methods. Figure HM shows that for p = 1,a@=0.1,b=1,c=0.5 
for x(t), relative errors between numerical schemes and exponential matrix 
method in [BY]. All figures denoted the results of the simulation based on 
equation systems (2), (3), (14), and (1) the values are used as mentioned 
above. Equation systems (2), (23), (14), and (4) correspond to equations 
in (Q). 

We compared CPU times of numerical methods on the Table I. Here, we 
used N = 1000 for all methods except Runge Kutta methods that we used 
N = 50. As we have seen in the Table [I], we can say, in case each explicit 
and implicit methods evaluated between oneself has not extremely difference. 
However we cannot say the same for the Runge Kutta. On the Table 2, we 
denoted qualitative results of the fixed point F2 for different time step sizes 
for N = 1000, p = 1, (20, yo) = (25,15), and K = 40. 
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26 19 26 19 
xt) 
24 yy | 418 24 18 
22 17 22 7 
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18 15 18 15 
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1" 10 —_—___— 1" 
10 8 10 
15 20 25 ° 5 10 15 20 25 


x(t) s x(t) - 
y(t) 18 y(t) 18 
(c) (d) 
Figure 1: For p=1, Theta and NSFD Schemes in explicit and implicit form: 
(a) Theta-explicit (b) Theta-implicit (c) NSFD-explicit (d) NSFD-implicit 
30 = 30 = 
yo yt) 
20 15 zo 15 


30 20 30 20 


x) 
vit) 


20 15 20 15 


a 
| 


Figure 2: For p=0.6, Theta and NSFD Schemes in explicit and implicit form: 
(a) Theta-explicit (b) Theta-implicit (c) NSFD-explicit (d) NSFD-implicit 
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xt) 
yo} 418.5 


xi) 
vit) 


26 19 25 18 
«) x(t) 
24 yo] jiss 70) 
18 
20 7 
175 
7 
16.5 
15 16 
16 
<— 7 et aes 
15 10 15 
25 ° 5 10 15 20 25 


Figure 3: For p=0.3, Theta and NSFD Schemes in explicit and implicit form: 
(a) Theta-explicit (b) Theta-implicit (c) NSFD-explicit (d) NSFD-implicit 


14 14 
13 13 
12 12 
"W " 
10 10 
8 10 12 14 5 20 22 24 26 8 10 12 14 8 20 22 24 26 
x(t) x(t) 


14 14 
13 13 
12 12 
"1 " 
10 10 
8 10 12 14 5 20 22 24 26 8 10 12 14 8 20 22 24 26 
x(t) x(t) 


Figure 4: Phase-portraits for p = 1, N = 2000, a 0.1, 6 Lac 0.5, A 0.01, 
(xo; Yo) oe (25, 15), k = 40 
(a) Theta-explicit (b) Theta-implicit (c) NSFD-explicit (d) NSFD-implicit 
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0.18 Be 
046 \ 


T ——r 7 T . . T 
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Figure 5: For p=1,a=0.1,b=1,c¢=0.5, (#o0,yo) = (10,10), K = 20, for x(t), 
relative errors between numerical schemes and exponential matrix method in [BY] 

— — —— EMM and Theta-explicit 

EMM and NSFD implicit 

sacceeeecaane EMM and Theta-implicit 

ao. EMM and NSFD-explicit 


Table 1: CPU times (seconds) for (ao, yo) = (25,15), h = 0.01 


p Theta Theta NSFD NSFD Runge Kutta 
Explicit | Implicit | Explicit | Implicit 
0.3 | 1.912602 | 2.216342 | 2.259337 | 1.813677 3.597516 
0.6 | 1.909790 | 2.522203 | 2.331288 | 1.77084 3.6738 
1 1.86862 | 2.273820 | 4.053474 | 1.767188 4.222919 


Table 2: Qualitative results of the fixed point E for different time step sizes for 
N = 1000, P = 1,(x0, yo) = (25,15), kK = 40(Con. and Div., resp., Convergence and 


Divergence) 
h Theta Theta NSFD NSFD | Runge Kutta 
Explicit | Implicit | Explicit | Implicit 

0.01 Con. Con. Con. Con. Con. 

0.1 Con. Con. Con. Con. Con. 
1 Div. Con. Div. Con. Div. 

10 Div. Con. Div. Con. Div. 

100 Div. Con. Div. Con. Div. 
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7 Conclusion 


In this paper, fractional-order form of the Hantavirus model is introduced. 
The stability of the equilibrium points is studied. Two numerical methods 
have been presented in the explicit and implicit form for solving fractional 
order Hantavirus epidemic model. When we use implicit solutions, we need 
to use the Newton—Raphson method to solve implicit system by converting 
to explicit form. For NSFD schemes, one denominator function and different 
nonlocal terms have been proposed, and the results have been compared 
with each other. As it seen clearly that explicit and implicit NSFD methods 
should be used ultimately depending on the choices of the nonlocal terms. 
On the other hand, explicit methods produce the same accuracy, but with 
less computational effort and time than implicit methods. 
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